Magnetically amplified light-shining-through-walls via virtual minicharged particles 



m 

o 

(N 

m 



Oh. 
I , 

CD 



(N 
> 

00 
On 
^, 

cn 

o 

(N 



X 




Babcttc Dobrich,^'^'G Holger Gics,^'^ Norman Ncitz/'[j and Felix Karbstcin^'^ 

^ Theoretisch-Physikalisches Institut, Friedrich-S chiller- Universitdt Jena, Max-Wien-Platz 1, D-07743 Jena, Germany 

^ Helmholtz-Institut Jena, Frobelstieg 3, D-07743 Jena, Germany 

(Dated: January 24, 2013) 

We show that magnetic fields have the potential to significantly enhance a recently proposed 
light-shining-through- walls scenario in quantum-field theories with photons coupling to minicharged 
particles. Suggesting a dedicated laboratory experiment, we demonstrate that this particular tun- 
neling scenario could provide access to a parameter regime competitive with the currently best direct 
laboratory limits on minicharged fermions below the meV regime. With present day technology, 
such an experiment has the potential to even overcome the best model-independent cosmological 
bounds on minicharged fermions with masses below C'(10~*)eV. 

PACS numbers: 14.80.-j, 12.20.Fv 



I. INTRODUCTION 

"Light-sliining-through-wall" (LSW) experiments have 
become a valuable tool in the search for numerous the- 
oretically well-motivated 'weakly interacting slim par- 
ticles' (WISPs) that could exist beyond the standard 
model of particle physics. Here, the shorthand 'WISP' 
subsumes hypothetical light particles at the (sub-)eV 
scale which exhibit a weak effective coupling to the 
electromagnetic field and are thus sensitive to optical 
searches [l|, 01 . Such particles generically arise in a num- 
ber of extensions of the standard model and under certain 
conditions also constitute a viable dark matter candi- 
date 0^9 • Most prominently, LSW experiments provide 
for the best direct laboratory bounds on axions, axion- 
like particles (ALPs) jTHlSj and massive hidden-sector 
photons (also referred to as paraphotons) [Tj], as well 
as for indirect limits on massless hidden-sector photons 
and fractionally charged or 'minicharged' particles |15l - 
Il7| . See, e.g., [l^l and references therein for a recent 
overview. 

In classic LSW scenarios the barrier is 'tunneled' by 
real, i.e., on-shell WISPs. Their paradigm is the LSW 
scenario with ALPs [1^, [201 ■ Light-shining-through- walls 
becomes possible if the probe photons are converted into 
ALPs in front of the wall and are reconverted into pho- 
tons behind the wall. In typical laboratory searches the 
conversion processes are induced by strong dipole mag- 
nets, see, e.g., [1, [ifl [ill [l3, [2l| ■ We depict this scenario 
in Fig. [H 

Besides LSW experiments, polarimctric measurements 
[22, 23] carried out in various setups [24l426| constitute an 
important means in the general search for WISPs. Re- 
cent suggestions involve, e.g., advanced interferometric 
setups [27H29J . tests of Coulomb's la w jBOl. 311 or experi- 



ments utilizing high-intensity lasers [32|-[34 



* Now at DESY, Notkestrafie 85, D-22607 Hamburg, Germany. 
t Now at Max-Planck-Institut fiir Kernphysik, Saupfercheckweg 1, 
D-69117 Heidelberg, Germany. 



FIG. 1. Classic LSW scenario with ALPs. In the presence of 
an external magnetic field (represented by wiggly lines ending 
at crosses) a photon entering from the left hand side (wiggly 
line) can be converted into an ALP (dashed line). Due to 
its weak coupling to matter, the ALP can traverse a light- 
blocking barrier nearly unhindered and can subsequently be 
reconverted into a photon. 



In addition, also LSW experiments using astrophysi- 
cal sources, particularly the sun, have been carried out 
with remarkable success as, e.g., in the CAST exper- 
iment [33, see also [Sa]. Such LSW 'helioscopes' as 
well as arguments involving stellar energy loss observa- 
tions, lead to extremely strong constraints [37H39J . The 
drawback is that most of these constraints are somewhat 
model dependent, as their derivation involves extrapola- 
tions to a different momentum-transfer regime [401 [4l| . 
Further indirect bounds on various hypothetical parti- 
cles can be derived from cosmology by studying possible 
distorsions of the cosmic microwave background (CMB) 
[43 |. The model- independent limits of these bounds for 
minicharged particles are currently somewhat stronger 
than the direct laboratory bounds. 

The high sensitivity of laboratory LSW experiments 
stems from the large number of laser photons available 
for the conversion process in front of the wall, facing the 
possibility of even single photon detection on its rear 
side. Currently, the best experimental bounds are due 
to the first stage of the 'any-light particle search' ALPS- 
I [1^ [43 at DESY. Proposals to enhance the sensitivity 
of this experiment even further, involve an enlarged mag- 
netic length |4J] and the installation of a second cavity 
on the regeneration side, i.e., utilizing so-called 'resonant- 




FIG. 2. LSW scenario with virtual minicharged particles, 
also referred to as 'tunneling of the 3rd kind', cf. [52]- A 
photon can traverse a light-blocking barrier through a virtual 
fermionic or scalar particle-antiparticle loop. This process is 
possible at zero field, cf. [52|. Here we study the process 
in the presence of an external magnetic field. The dressed 
minicharged particle propagator, involving an arbitrary num- 
ber of external field insertions, is represented by the solid 
double line. 



regeneration' |45l - l48| . Higher energy photons [12|, |49| are 
used in other contexts. Noteworthy, the latter concept 
could recently be verified experimentally even for regen- 
eration in the sub-quantum regime [50| . 

In this paper we aim at investigating a different LSW 
scenario, namely the 'tunneling' of a barrier through in 
general virtual particle-antiparticle intermediate states. 
Our main focus is on the process in the presence of an 
external magnetic field, cf. Fi g. [^ A compact version of 
our results has been given in |51| . 

As the tunneling process with the virtual particle- 
antiparticle pair complements both quantum mechanical 
tunneling and the tunneling of real particles as depicted 
in Fig.[l] it has been baptized 'tunneling of the 3rd kind' 
[53 | . Although such a process is in principle also possible 
with neutrinos, this standard model background is highly 
suppressed due to the Fermi scale. 

Although the photon-to-photon transition amplitude 
grows logarithmically with decreasing mass of the 
minicharged particles in the absence of an external field, 
zero-field 'tunneling of the 3rd kind' cannot improve ex- 
isting bounds of the minicharged particles' parameter 
space [53 . By contrast, we show in this work that current 
laboratory bounds on minicharged fermions in the sub- 
meV regime can be significantly enhanced by an external 
magnetic field. Most importantly, they can even out- 
match the best current model-independent cosmological 
bounds below minicharged fermion masses of C'(10^"')eV. 

Our paper is organized as follows: In Sect. |ll]the gen- 
eral formalism describing the LSW scenario depicted in 
Fig. [2] is outlined and the photon-to-photon transition 
probability is derived. Sect. IIIII focuses on the explicit 
determination of the transition probability. This in par- 
ticular requires an adequate analytical treatment of the 
photon polarization tensor. Appropriate limiting cases 
which ultimately allow for a numerical evaluation of the 
transition probability are worked out. In addition, an- 
alytic asymptotics in the most relevant limits are de- 
rived. Specializing to a specific, experimentally feasible 
setup, in Sect. IIVI the corresponding exclusion plot in 
the fractional-charge mass plane are worked out and dis- 
cussed in detail. The paper ends with conclusions and 



an outlook in Sect. M 



II. SETTING 

Let us briefly review the basic setting as introduced 
in [52J | and extend it to account for an external magnetic 
field B 7^ 0. We adopt the metric convention 5^1, = 
diag(— 1, +1, +1, +1) of [52[, such that the momentum 
four- vector squared reads k^ = k^ — w^. 

We start with the effective Lagrangian describing pho- 
ton propagation in a constant external magnetic field of 
strength _B = |B| given by. 



1 



C[A] = --F,^ix)F^-^ix) 



A^{x)n'^''{x,x'\B)A,{x') , (1) 



where 11'*'' (x, x'\B) denotes the photon polarization ten- 
sor in the presence of the external field, F^i^ the field 
strength tensor of the propagating photon A^, and x a 
spatio-temporal four-vector. 

As a result of their negligible coupling to matter, the 
wall does not affect the propagation of minicharged par- 
ticles. Hence, translational invariance is maintained for 
the minicharged particles and thus for the polarization 
tensor at one-loop order. As long as the magnetic field is 
homogeneous in the relevant spacetime region, this im- 
plies nf"'{x,x'\B) = Il>"'{x-x'\B). 

Upon a variation of Eq. ([1]) and a transformation to 
momentum space, the following equation of motion is 
obtained 

(/c V - kf'k" + Il^'''{k\B)) A^{k) = . (2) 

The photon-to-photon transition probability depends on 
the polarization mode of the photons. Contrary to the 
zero-field situation, there are three independent polariza- 
tion modes in the presence of an external field, which we 
label by an index p = 1,2,3. As the vacuum speed of 
light in external fields deviates from its zero-field value 
[53I455J ■ and the vacuum exhibits medium-like properties, 
the occurrence of three (instead of two in the absence of 
an external field) independent polarization modes is not 
surprising. For instance, the corresponding polarization 
effects induced by minicharged particles have been stud- 
ied in [13, [5^. Further diffractive effects in inhomoge- 
neous fields as also considered in quantum electrodynam- 
ics (QED) [531 are conceivable and certainly worthwhile 
to be studied. 

Let us now bricfiy introduce the explicit structure of 
the photon polarization tensor in the presence of an ex- 
ternal field. Here we limit ourselves to the special case 
of a purely magnetic field [S^ . This naturally suggests a 
decomposition of the photon four-momentum fc^ in com- 
ponents parallel and perpendicular to the magnetic field 
vector B. The angle <);(B,k) is denoted by 6. Without 



loss of generality B is assumed to point in ei direction, 
and the following deeomposition is adopted, 



nJ — Ku |~ 1^ I ; 



iO,0,k^,k^). 



(3) 



In the same manner, tensors can be decomposed, e.g., 
9^'^ = 9u'^ + qT" ■ It is then convenient to introduce the 
following projection operators onto photon polarization 
modes, 



pfj.v 

^1 



k'^K 



and 



pfj.1^ 

^2 



01 



ll 



Defining a third projector as follows. 



Ps'" = g'" 



kf'k" 

fc2 



pfj,u pfJ.l> 



(4) 



(5) 



the three projectors P^" obviously span the transverse 
subspace. Whereas two of these polarization modes can 
be continuously related to the photon polarization modes 
in the absence of an external magnetic field, the third 
mode manifests itself in the presence of the external field 
only. 

Note that P|"^ and P2"' have an intuitive interpreta- 
tion, provided that k jj" B. Namely, they project onto 
photon modes polarized parallel and perpendicular to the 
plane spanned by the two vectors, k and B. For k jf B 
these are the polarization modes, that can be continu- 
ously related to those in the zero-field limit. Remark- 
ably, for the special alignment of k || B, the situation is 
different. Here, the modes 2 and 3 can be continuously 
related to the two zero-field polarization modes. 

With the help of Eqs. dU and ([5]), the photon polariza- 
tion tensor in a purely magnetic field can be decomposed 
as follows [5i,[6a], 



n^''(fc)= ni(fc)Pf 



U,ik)Pr 



U,ik)Pr 



(6) 



where the scalar functions Up are the components of the 
polarization tensor in the respective subspaces. Thus, 
dropping Lorentz indices, the equation of motion for pho- 



tons in mode p, 
reads 



P A" 



resulting from Eq. ^ 



{k^+Ilp{k\B))Ap{k) = 



(7) 



As depicted in Fig. [S] we assume the probe photons to 
propagate in e^ direction, normal to the light blocking 
barrier, i.e., k || e^. The barrier extends from r = to 
r ^ d along e^, and has infinite extent orthogonal to e^. 

Introducing partial Fourier transforms of the photon 
field and the polarization tensor. 



Ap{r,cu)^ J^e^^'^Apik), 



(8) 



Up{r-r',io\B) 



dk 
2^ 



g,(r-r')kjjp(^|^) ^ (9) 




FIG. 3. The probe photons are assumed to enter in e^ di- 
rection normal to the wall. We look for regenerated photons 
behind the wall propagating along the line of propagation of 
the incident photons. The angle between k and B is denoted 
bye. 



where k = |k|, the equation of motion becomes 



{uj^ + d^)Ap{r,uj)^ dr'np{r~r',uj\B)Ap{r',Lo)., 



= :jp{r,uj\B) 

(10) 
where we have defined the fluctuation-induced current j. 
Following j52| , we choose reflecting boundary conditions 
on the front side of the barrier, r = 0, for the incident 
photons. Note that the calculation could easily be gen- 
eralized to other boundary conditions. Accordingly, the 
induced current within and beyond the barrier reads 



jpir > 0,tj\B) = / dr"np{r~r",iu\B)a{uj) sin(a;r") , 

(11) 
where a(a;) denotes the amplitude of the incident pho- 
tons. Moreover, the photon source is assumed to be suf- 
ficiently far away, such that the lower integration limit 
in Eq. (|lip can be formally sent to — (X) (cf. comment at 
the end of this section). 

To obtain the photon-to-photon transition probability, 
the outgoing photon wave on the rear side of the barrier 
has to be determined. For detector positions asymptoti- 
cally far beyond the barrier at r 3> d, we find 



'ir:^d,uj\B) 



CK3 

i f dr' 

d 



^iu>{r—r') 



2UJ 



■ jp{r' ,uj\B) , (12) 



where we formally sent the upper integration limit to 
+00, employed the free outgoing Green's function for 
the operator [uP 4- 9^) ^"^^ restricted ourselves to the 
right- moving, i.e., transmitted components of the outgo- 
ing photon field. This corresponds to absorbing bound- 
ary conditions on the rear side of the wall at r = d. 
Thus, combining Eq. pT|) and Eq. P^ . the photon-to- 



photon transition probability can finally be written as 



P,7-i-7 



lim 



Ai^^r.uj) 



a(uj) 



dr' ■ 



2uj 



dr"Up{r' ~r")sm{ujr") 



(13) 



where we have dropped any explicit reference to the mag- 
netic field B for clarity. 

As lip (r'—r") receives its main contributions from rel- 



ative distances \r'^ 



of the order of the Compton wave- 



length ~ 1/to of the minicharged particles, and falls off 
rapidly for |r' — r"| ^ 1/?ti, it is justified to formally ex- 
tend the integration ranges to ±cx), respectively. With 
respect to an actual experimental realization this implies 
that the magnetic field has to be sufficiently homoge- 
neous only within a sphere of diameter > 1/m centered 
at the intersection of the optical axis with the wall. 

In the subsequent section we explicitly evaluate 
Eq. (J13p . Our focus is on minicharged Dirac fermions. 
Remarks on the rather different physics of scalar 
minicharged particles are postponed to App. \^ 



III. POLARIZATION TENSOR AND 
TRANSITION PROBABILITIES 

A. Photon polarization tensor in momentum space 

In the presence of an external electromagnetic field, the 
photon polarization tensor in momentum space is explic- 
itly known at one-loop accuracy. Utilizing the proper- 
time representation [55|, it is most convenie ntly stated 
in terms of a double parameter integral [59|, |61[, which 
in general cannot be tackled analytically (see, e.g., [S^ 
and references therein). Here we are interested in the 
special case of a purely magnetic field [S^ only. Several 
well-established approximations to the polarization ten- 
sor in a purely magnetic field have been worked out [62l - 
|65|. However, most approximations considered so far put 
special attention to on-the-light-cone dynamics. Besides, 
they are in general limited to physical problems that can 
be treated directly in momentum space, as their deriva- 
tion involves constraints to a certain momentum regime. 
In contrast, our approach involves a Fourier transform to 
position space, cf. Eq. (j9]), and hence manifestly requires 
knowledge about the full momentum regime. 

To achieve this, we turn to the explicit expression 
of the photon polarization tensor in a purely magnetic 
field at one-loop accuracy (as reproduced in our nota- 
tion in [60| ) . The photon polarization tensor of standard 
quantum electrodynamics is adapted to minicharged par- 
ticles by substituting the electron-positron loop for a 
minicharged particle-antiparticle loop. Formally, this 
amounts to a change of the coupling e — >■ ee, with e 
referring to a dimensionless fractional coupling, and an 



identification of m with the mass of the minicharged 
particles. Noteworthy, higher loop corrections to the 
photon polarization tensor at one-loop accuracy adapted 
to minicharged particles are therewith suppressed by at 
least an additional factor of e^. 

We perform formal expansions in eeB — > and eeB — > 
oo. The combined quantity eeB appears as the natural 
parameter, since the coupling of the external field B to 
the minicharged particle line is ee. The corresponding 
(dimensionless) physical expansion parameters, and the 
ranges of applicability of these expansions will be dis- 
cussed on the level of the photon-to-photon transition 
amplitude below. 

A perturbative expansion of the photon polarization 
tensor in powers of eeB, while retaining its full momen- 
tum dependence, is straightforward. It results in 



+ 0{ieeB)^), (14) 



were the upper index (n) refers to the order in the expan- 
sion, i.e., denotes contributions of order (eei?)". In con- 
sequence of Furry's theorem, the expansion is in terms of 
even powers of eeB only. H^"^ is the photon polarization 
tensor in the zero- field limit. 
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Tl, 


= n(o) + 


n^^ 


Us J 




[< 



n(o) = (ky 



in 



diy 



V 
$0 



with 



<I>o = 771 — irj + 



l-u^ 



(15) 



(16) 



Here 77 > denotes an infinitesimal parameter and a ~ 
e^ /{Ait) is the fine-structure constant. The contribution 
at O {{eeBf) reads 



{1-uy 




(17) 



Carrying out the expansion eeB — >■ cxj, while keeping 
the full momentum dependence, turns out to be more in- 
volved. To leading order in the parameter eeB we obtain 



e^aeeB 




dv < 








where we have introduced [cf. Eq. (fTB]) ] 



$n 



11] 



(18) 



(19) 



Sublcading corrections to Eq. ((18)) arc at most logarith- 
mic in eeB. Hence, the strong field limit is dominated by 
ni(fc). 112 (fc) and n3(fc) are suppressed and start con- 
tributing at subleading order. Note that a factorization 
with respect to the momentum dependence, kn^ and /c^, 
is encountered in Eq. p^ . 

Equations, p^ and (|18p are the representations of the 
polarization tensor which serve as the basis of our further 
analysis. 



B. Towards the polarization tensor in position 
space 

Let us now explicitly implement the one-dimensional 
Fourier transform, Eq. ([9]), of the polarization tensor 
from momentum to position space. We first focus on the 
regime where the pcrturbative expansion in eeB given in 
Eq. ([13]) holds. The Fourier transform of Eq. dH]) in- 
volves terms of the general form 



dk 

2^' 



,?(r— r')k 



= 4' 




\ u;^ 1 — u-^ 



1^ (20) 



where -P(k) denotes a polynomial in k, and I E N. Equa- 
tion (|13p involves further integrations over the spatial 
coordinates in the argument of n(r — r') such that the 
difference r — r' is always non-zero and strictly positive. 
The evaluation of the k integral with the Residue the- 
orem requires special care, as the location of the poles in 
the complex k plane depends on the further integration 
parameter v. Noteworthy, for given to > 2m, the quan- 
tity uj^ ~ T^i;^ exhibits a sign-change as a function of i^. 
Defining |52| 



A= Wl- 



4m2 



a;2(i_i.2) ' 



(21) 



and 



Am? 



a;2(l-i.2) 



(22) 



the poles of the integrand in Eq. ((20)) in the complex k 
plane are located at 



k = ±iwK ± (1 -f 1)77 for oJ^ < 

k = ±a;A ± (1 + i)-q for o;^ > 



4to2 

47^2 



(23) 



This suggests the following decomposition of the v inte- 
gral in Eq. ((T4| . 



dv h[v) 



^4" 



A-^ 


^2 
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dA 


I"'"' 







00 
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dti 


-T-h{v 
dn 



^(l+»^) 



(24) 



where h{v) denotes the integrand of the v integral. With 
these preparations, the Fourier transform of Eq. ((H)) is 
implemented straightforwardly. Subsequently we rely on 
the decomposition in Eq. ((^^ . It is possible to distin- 
guish two cases: For a; < 2m, there is only the k integral. 
For cj > 2?7i, both the k and A integrals contribute. By 
virtue of the pole structure in the momentum integration, 
cf. Eq. ([23]) J the integrand of the k part is exponentially 
damped with increasing k, whereas the A part is oscillat- 
ing as a function of A. Henceforth we label contributions 
due to K by "<", and contributions due to A by ">". 

In the strong field limit with large eei?, Eq. ((T5)) . the 
Fourier transform involves terms of the form 



dk^,(,_,')kl^^.2 -, 



27r 



*0||| " 

dk e*('-'^')k (k2 



2 jj.2 gj^2 Q 



2tt 



k2- 



cos-^ 6 \ uj'^ 1 — iy'^ 



(25) 



Due to the factor quadratic in k in the exponential, the 
integrand does no longer vanish in the limit k — >■ ioo, 
and the evaluation of the k integral becomes more com- 
plicated. It is most conveniently performed with the help 
of the convolution theorem, 



dk 

2^ 



c^'^'-gi(k).g2(k)EE / dr~.9i(f)g2(r-f), (26) 



where 



Mr) 



CX3 

y'ge*^g.(k) (* = 1,2). (27) 



entifying 




k2 "^ 


(28) 


V2 <^^ ( 1 4m2 1-11, \ 


k^ sin2 e 


(29) 



the Fourier transform of gi (k) resembles that in Eq. ()20p . 
It can be performed with the Residue theorem. The de- 
composition of the V integral, Eq. (|24l) . remains favorable. 



but the poles in the complex k plane are shifted to 

4m^ 



^'^i^^^'^'^^ 


for 


0J'< 


k-± "^^±(1 + 1)^ 


for 


w2> 



COS0 



4771^ 
l-:/2 



(30) 



The function 52 (k) does not depend on v. Its Fourier 
transform reads 



.92 (r) 



eeB 



2-K sin^ B 



(31) 



Therewith, the corresponding convolution integral on the 
right-hand side of Eq. (|26p can be rewritten in terms 
of error functions. Let us however note that Eq. ()3ip 
approaches the Dirac delta function in the formal limit 
^iS-^ — > 0, via the heat kernel 



,,2 



<5(r) 



(32) 



This implies, that in the strong field limit, we are inter- 
ested in here, we may approximate g2(f') — '^(t'). Using 
this in Eq. ((26)) . the Fourier transform to be implemented 
in Eq. ^^ reduces to that of Eq. (f28| . 



C. Tow^ards the photon transition amplitude 

Next, we perform the two position-space integrals in 
Eq. p^ . This can be achieved with the identity 

c« 

dr' I Ar" [i(/-r")]'e'^^^('''-''"^sin(wr")e-'"'^' 

d —00 



2+/ 



(^a)' 



,(_A^_l)d 



(cos B-iT]- A)2 (cos 6 + X) 



(33) 



where / e Nq. Note that we explicitly account for the 
irj prescription only in such terms where it is essential 
for the further calculation. The corresponding equation 
for K follows by substituting A — > in. After this sub- 
stitution, the poles induced in the complex k plane are 
purely imaginary. For angles < ^ < f , not too close to 
9 = ^, these poles are located reasonably far away from 
the interval of integration, which corresponds to the real 
K axis. In case of the A contribution, the situation is dif- 
ferent. Here Eq. ([55)) gives rise to poles at A = —cos 61 
and A = cos 6 — irj. If the contour of the A integration is 
chosen along the real A axis, and 9 is restricted as above, 
the first pole always lies well outside the interval of inte- 
gration. However, the latter pole touches the integration 
contour if 



cos6l<5R\/l 



4,m? 



O uj sin 9 > 2m . 



(34) 



Let us remark that, sticking to a perturbative expan- 
sion in powers of eeB, and assuming non- vanishing — , 
Eq. ([33]) is employed for = only. This is a consequence 
of the fact that all the Residues in Eq. p^ arise from 
inverse powers of Eq. (fT6)) . As condition (p4| cannot be 
fulfilled for 9 = Q and — 7^ 0. it is of no relevance in the 
perturbative regime. 

It is convenient to introduce dimensionless auxiliary 
functions /< and /> as in [52[, and to decompose the 
photon-to-photon transition probability in Eq. (J13p ac- 
cordingly. 



p , — 



367r2 



./<+/> 



(35) 



The modulus in Eq. (f55]) allows us to drop global phase 
factors in the definition of the auxiliary functions in the 
remainder. 

Before turning to the situation of a non- vanishing ex- 
ternal magnetic field, let us recall the auxiliary functions 
at zero field. As already derived in [52l |. for transversal 
photons we find 

„ — U>dK 



/f = 



K^ 



dK ■ 



i + K 



,.2 4m^ 
rL — T- 



(1 + ^2)3/2 



(36) 



/f = 



1^1 



„iujdX 



dX- 



l-A 



1-A2 



4m^ 



1-A2 + 



2nr 



(37) 



(1 - A2)3/2 

The corresponding transition probability at zero field 
reads 

(38) 



p(0) = li^ 



4°' 



/^°' 



In addition to the results of [52l | , we can provide for ana- 
lytic asymptotics of Eq. (|38p in the limit where a; ^ 2m 
and d < ^i;- {'^y^ (cf. App. E]). Starting from A == 0, 
the integrand in Eq. (J37]) rapidly increases with A. The 

factor A /l — A2 — ^^ guarantees that it is ultimately 
bent back and vanishes exactly at the upper integration 
limit. This gives rise to a peak close to the upper inte- 
gration limit. The smaller the ratio ^^, the closer the 
upper integration limit approaches the pole of the inte- 
grand at A = 1, and the more pronounced is the peak. 
This motivates to approximate the following part of the 
integrand as. 



iiudX 



(l_A2 + 22|l)e 

(l-A)(l-A2)3/2 

and the square root by 



1- A+^ 
V2(l - A)5/2 



1-A2- 



V2 



-A 



1/2 



(39) 



(40) 



This approximation retains the important feature that 
the integrand vanishes at the upper integration hmit, a 
behavior that would inevitably be spoiled by a naive se- 
ries expansion of the integrand in Eq. I\'37\i about A = 1 . 
The resulting expression can be integrated explicitly. 
Combining Eqs. ([39| and (|40|). we obtain 
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(41) 



The integrand of the k contribution, Eq. ([55)1 , does not 
feature any pole close to the interval of integration and 
is highly suppressed as compared to Eq. ([37|) . Moreover, 
it decays exponentially with cudn. Hence, in the limit 
^ < 1 and for d < ;ij (ir)"^ ("^f- App. E]), we obtain 



/r = -21n('^ 



^iujd 
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' 2m •] 



f^' 



o (tt) 



(42) 



and Eq. 



is well approximated by 



p(zcro)^f^l^2/2mA 
^7^7 - 9,^2 '"^ {oj J 



(43) 



which is in agreement with the numerically determined 
asymptotics presented in Eq. (3.17) of j53 |. 



D. Transition amplitude for •weak magnetic fields 

We now turn to the situation of a finite magnetic field. 
It is instructive to determine the photon-to-photon tran- 
sition probability in the pcrturbative regime first. The 
leading perturbative correction to the zero-field limit is 
depicted schematically in Fig. 01 

In this section, we limit ourselves to an angle of 6* = 0. 
This significantly simplifies the expressions of the respec- 
tive auxiliary functions, without neglecting any of their 
basic features. The results for non-zero angles follow 
straightforwardly from Eq. p7)) . 

As elaborated in Sect. |lTl in this situation the modes 
labeled by p = 2, 3 correspond to the two photon prop- 
agation modes that can be continuously related to zero- 
field polarization modes. Due to rotational invariance 
along the magnetic field vector, the components of the 
polarization tensor for these modes coincide for 9 = 
[cf. Eq. pT)) ]. With the preparations in Sects. HITB] and 
nil C[ it is straightforward to evaluate the corresponding 



transition amplitudes. At order (eeB)"^, one finds 

2 
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W^ 



V^ 
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i + 2k. + dujK{i + k) 



K^{i + k)2(1 + K2)3/2y'l + K2 



4?-n^ 



(44) 



The corresponding function /> can be obtained from /< 
by substituting k — > —iX and changing the integration 
boundaries accordingly [cf. Eq. p4))]. This yields 



./: 



(2)^ /eeSy 4m2 
\uj^ ) a;2 



K^I 



dAe*"''^ 



1 - 2A - ida;A(l - A) 



A2(l - A)2(l - A2)3/2y'i_A2- 



(45) 



We briefly comment on the ranges of validity of these ex- 
pressions. A perturbative expansion in powers of eeB is 
trustworthy when the dimensionful quantity eeB is small 
in comparison to all the other quantities of the same di- 
mension available. For the setting considered here, this 
means that the following dimensionless ratios involving 
eeB are small: 



eeB 



< 1 and 



eeB 



<1 



w^ 



(46) 



Note that we do not include a ratio involving the length 
scale d in Eq. (|46)) . This is due to the fact that d never 
appears alone, but always in the dimensionless combina- 
tion ujd. 

Let us have a closer look at Eqs. (j44]) and (|45p . In 
contrast to the zero-field limit, where Eq. ((36|) is always 
finite, Eq. pi)) is found to diverge for to = 2m and all val- 
ues of d. As this divergence or 'resonance phenomenon' 
is absent in the zero-field limit, it can be considered as a 
genuine manifestation of the external field. Note that 
it has already been discussed by previous authors on 
the level of the polarization tensor of QED, see, e.g., 
[64, 68 70] . However, it becomes particularly relevant 
for minicharged particles, since, in contrast to QED, the 
size of m could easily be smaller than the photon energies 
employed. 

The divergence in the photon-to-photon transition 
probability for w = 2m of course signals a break down 
of unitarity in our calculation, as it would predict an ar- 
bitrarily large number of outgoing photons for any small 
number of incoming photons. This unitarity violation 
is a consequence of the idealized limit of a perfectly co- 
herent infinite incoming wave train. It would require a 
proper treatment, e.g., by taking the finite line width of 
the laser into account. Yet, as the resonance may indeed 
lead to a strong sensitivity to minicharged degrees of free- 
dom fulfilling the condition lo = 2m, a careful analysis 
seems highly worthwhile. However, in the present work 
we continue to work in the perfectly coherent wave limit 



and ignore potential enhancements arising from such res- 
onances for our phenomenological conclusions, i.e., focus 
on the parameter space that can be firmly excluded even 
if the resonances would be smoothed out in an actual 
experimental realization. 

A direct numerical evaluation of Eqs. (03]) and (pS)) 
becomes tedious for u > 2m. In this case the integrand 
in Eq. (|44p seems to diverge at the lower integration 
limit and that in Eq. (|45)) at both the lower and the up- 
per integration limits. However, the ir] prescription (cf. 
Sect. IIIIBp . assumed implicitly in Eqs. (|44]) and (|431) . 
ensures that these divergences do not lie on the integra- 
tion contour. Apart from the true divergence at w = 2m 
discussed above, the encountered superficial divergences 
can be circumnavigated with the help of integration by 
parts. 

Using the following identity for cj ^ 2m, 
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1 + k2 



^t^/l 



4m2 " 



'C, (47) 



with C denoting an integration constant, Eq. (^l)) can be 
rewritten as 
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4m-' 



(2) ^ f eeB 

i^+ujd[ujd{l + K'^)^2i + 
(i + «;)(1 + k2)5/2 



dK ( 



5k 



(48) 



and a surface term. Note that the lower label in Eq. 
is "<" , rather than "<" , as the case w = 2m has been ex- 
plicitly excluded. Moreover, we have omitted the surface 
term of the integration by parts in Eq. (|^5)) . It vanishes 
identically for oj < 2m. For ui > 2m, it cancels with the 
corresponding surface term arising in an analogous inte- 
gration by parts of Eq. (|^5]) . This can be easily verified 
by taking into account the explicit irj prescription (cf. 

Sect.MM- 

Correspondingly, employing Eq. (j47p with n -> — iA, 
and skipping the surface term, Eq. ps)) becomes 



(2) _ feeB 
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(49) 



Equations (^5)) and ([^^ are free of divergences in the 
intervals of integration, and are thus perfectly suited for 
a numerical evaluation. 

Finally, we combine the zero-field contributions, Eqs. 
([55)1 and (|37p . with the leading perturbative corrections. 




FIG. 4. LSW scenario in the perturbative regime with two 
couplings to the external field. This process merely serves 
as an illustration: as shown in Sect. IIIIDI the perturbative 
approximation ceases to hold far before the perturbative cor- 
rection adds significant contributions to the zero-field result. 



Eqs. P5|) and (|49l) . to obtain the photon-to-photon tran- 
sition probability, Eq. (j35p . in the presence of weak mag- 
netic fields. 
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Following the same reasoning as in the absence of an 
external field, we can provide for analytic asymptotics of 
Eq. ([nni) in the limit where a; > 2m and ^ < ;^ (^)~ 
(cf. App.|B]). Now our focus is on Eq. P5|) . We approxi- 
mate the following part of the integrand as. 



f^+i^d [ujdil - A2) - i{2 + 5A)] 



udX 



(l-A)(l-A2)5/2 

and use Eq. ([jn]) for the corresponding square root ex- 
pression in Eq. ((49]). To leading-order in — ^ 1, we 
obtain 
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0{(^' 



(52) 



Hence, in the limit ^ < 1, ^ < 1 and d < i (^) \ 
Eq. (0^ is well approximated by 



f(2) _ 



2 /eeB 
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As in the zero-field case, Eq. (|48|) is suppressed in com- 
parison to Eq. (|53|) in the considered limit. It is already 
clear at this point, that in the regime of validity of the 
weak field approximation, the transition probability for 
u! ^ 2m is still dominated by Eq. (|42]) . which grows log- 
arithmically with — ^- 0. Therefore, the leading contri- 
bution to Eq. (1501) agrees with the leading contribution in 
the absence of an external field, and is given by Eq. (|43l) . 



E. Transition amplitude for strong magnetic fields 

Let us now aim at insights into the opposite regime, 
i.e., the hmit of strong magnetic fields. It is immediately 
obvious from Eq. ^TE\\ , that here the leading contribution 
arises from the p ~ 1 component of the photon polariza- 
tion tensor, whereas the components p = 2, 3 contribute 
at subleading order only. 

Due to the overall exponential damping factor '^ 
exp [-kl/{2eB)] in Eq. ^ (cf. also ll M, III), 
the maximum photon-to-photon transition probability 
attainable, is to be expected for small angles 9. Recall 
however, that in the strict limit 6 = the mode p = 1 
can no longer be continuously related to any of the pho- 
ton polarization modes in the vacuum (cf. Sect.|lT|. It is 
not clear at all, if and how this mode can be excited for 
9 — 0, cf., e.g., the controversial viewpoints in |72l. l73|. 
Hence, our main focus is on a preferably small, but fi- 
nite angle 9^0. As discussed in detail in Sect. |lll here 
the p = 1 mode corresponds to transversal photons with 
polarization vector in the (k, B) plane. It resembles an 
Alfven-like transversal mode. The strict limit 9 = can 
be seen as constituting the asymptotics of an ideal exper- 
iment, allowing for an arbitrarily small but nonvanishing 
angle. In contrast, the second (magneto-acoustic) polar- 
ization mode that can be continuously related to a pho- 
ton polarization mode at zero field receives no dominant 
magnetic enhancement. 

Recall, that in order to obtain Eq. p^ . we have con- 
sidered the formal limit eeB — )■ oo. We now specify the 
physical parameter regime, where the corresponding ex- 
pression for the photon-to-photon transition probability 
is trustworthy. This turns out to be the case, if the fol- 
lowing dimensionless ratios involving eeB are large [66|: 



eeB 



> 1 and 



eeB 



up' sin ( 



>1 



(54) 



As in Eq. (|46|. we do not include a ratio involving the 
length scale d in Eq. ([M|) . The length d never occurs 
alone, but always in the dimensionless combination ujd. 
In contrast to Eq. (pS)) . the second ratio in Eq. ([M)) in- 
volves a factor of sin 9 in the denominator. This is a 
direct consequence of the factorization with respect to 
the momentum dependence, encountered for the photon 
polarization tensor in the strong-field limit [cf. Eq. (fT8|) ]. 
Note that the conditions in Eq. ((54)) also have an intu- 
itive interpretation. A magnetic field results in Landau 
level quantization of the momentum components perpen- 
dicular to the magnetic field vector B. On the level of the 
photon polarization tensor, this implies a factorization of 
the loop integral into a sum over discrete _L momentum 
components, and a continuous momentum integration for 
the II components. The discrete _L momentum contribu- 
tions can be reinterpreted in terms of an effective, level 
dependent mass. Contributions of two successive Landau 
levels can be seen as differing in their effective masses 
squared by Atti^ r-^ eeB. 



Hence, the first condition in Eq. ([M]) ensures that the 
effective mass of particles in higher Landau levels is sig- 
nificantly larger as compared to the zeroth one, where 
the mass is unscreened. The latter condition in Eq. ([5^ 
ensures that the probability of real pair creation at higher 
Landau levels is small. 

As already discussed in Sec. IIIIBl for 2l|^ _> Q the 
Fourier transform of Eq. ([29| becomes a representation of 
the Dirac delta function. This limit is fully compatible 
with the parameter regime specified in Eq. (j54p . For 
J^fg ^ 1 an approximation of Eq. ([31]) by the Dirac 
delta function is well justified, and further simplifies the 
calculation. 



Following the steps outlined in Sects. IIIIBl and IIII CI 
for the p = 1 mode in the strong-field limit, we obtain 
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and 
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Let us emphasize that the structure of Eqs. (|55p and 
(|56p is different from the analogous expressions in the 
weak-field regime. Most strikingly, the factorization with 
respect to the momentum dependence in Eq. (|T8)) induces 
COS0 terms in the denominators [cf. also Eq. ([5(I|) ]. As 
discussed in Sec. IIII Cl a pole in the interval of integration 
is induced in Eq. (fSEj) . if the condition (|M)) is met. 

It is in general convenient to rewrite Eq. (|56p by em- 
ploying the same integration by parts as in the weak field 
regime, cf. Eq. (gT)). This yields 
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For the A contribution, this particular partial integration 
is helpful only if the integrand does not feature any poles 



in the interval of integration, i.e., if cos9 > 3?-i/l 



4m2 
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It results in 



and (1621). wc obtain 
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As in Scct. lIIIDl these integrations by parts arc only pos- 
sible for bj 7^ 2m, and the transition probability diverges 
for to = 2m [GJ, l68l - [7(]| | . The surface terms in both the k 
and the A contributions vanish by themselves. 

For the situation cos 9 < 5R\/l — ^^. it is more con- 

venient to rewrite Eq. (|55p in such a way that no explicit 
reference to the infinitesimal parameter 77 is necessary. As 
the respective expression is rather lengthy, we postpone 
it to App. O 

The photon-to-photon transition probability in the 
strong-field limit is given by 



p(strong 
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367r2 
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In accordance with the previous sections, we aim at 
analytical insights into the limit where cj ^ 2m and 



d«k{^y'- 



We first consider the regime, where 
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Hence, in the strong field limit |^^, — #^^t^| ^ 1, and 

for 2Hi ^ 1 and d < ;ij (tt)"' (cf- App.©, Eq. ([Ml) is 
well approximated by 
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Exactly as in the situation of vanishing and weak fields, 
Eq. (I57|) (the K branch) is suppressed in comparison to 
Eq. (|64|) (the A branch) in the considered limit. There- 
fore, the asymptotics for the photon-to-photon transition 
probability in the regime (|60| is independent of 9, and 
reads 
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It is also possible to extract an analytic asymptotics for 

uj ^ 2m and d ^ -^ (^) i'^ ^-'^^ regime complemen- 
tary to Eq. ([M)) . characterized by Eq. ([M)). This neces- 
sitates a careful analysis of Eq. (JC1[) , applying the same 
strategies as employed when deriving the other asymp- 
totics in the limit a; ^ 2m above. For the same reasons 
as discussed before, the k contribution Eq. (|57p is negli- 
gible. We find (cf. App. [C]) 



If we want to obey Eq. ([60]) together with uj ^ 2m, the 
angle 9 has to be very small. This motivates an expansion 
of Eq. (|58| in powers of 9, 
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As the integrand does not feature poles in the interval of 
integration, any explicit reference to rj has been omitted 
here. We approximate the following part of the integrand 
in Eq. dH]) by 
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4v^ (1 - A)V2 
and the square root as in Eq. (PUJ) . Combining Eqs. PO)) 
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and for small angles 9, 
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(68) 



Let us provide for an intuitive explanation of the con- 
dition wsin^ = 2m, separating the two regimes, defined 
by Eqs. (|34|) and ()60p . It is the on-set condition for real 
pair creation. To see this, it is illustrative to perform a 
Lorentz transformation along the direction of the mag- 
netic field vector, such that the parallel momentum com- 
ponent k|[ of the photons becomes zero, i.e., kn = 0. Vec- 
tor components orthogonal to the direction of the Lorentz 



11 



boost remain unaltered. In this particular reference sys- 
tem (denoted by '), the particles are still subject to a 
homogeneous external magnetic field. This comes about 
as a Lorentz boost in the direction of the magnetic field 
does not induce an electric field. However, the light-cone 
condition for photons now reads 



(fc'") =0, 



(69) 



and the on-set condition for real pair creation becomes 

fc'° = |k_L|=2m ^ u}sme = 2m. (70) 

Going gradually from the regime uj sin 9 < 2m beyond 
this limit, i.e., to uj sin 6 > 2m, we expect the photon- 
to-photon transition probability to drop, as real particle- 
antiparticle pairs are unlikely to reconvert into photons. 

Note that the photon-to-photon transition probability 
diverges at the real pair creation threshold uj sin 9 = 2m. 
An analogous unitarity violation has been encountered 
already for u = 2m. As discussed in detail in Sec. lIIIE)! 
it can be cured by taking into account the finite line 
width of the laser. Let us remark, that the divergence at 
a; sin 6 = 2m is not visible in any pcrturbativc expansion 
in powers of eeB. Perturbative calculations account for 
a finite number of external field insertions in the loop 
diagrams only, and the intermediate minicharged parti- 
cle lines are described by zero field propagators. Hence, 
translational invariance is not broken in the propagators 
[cf. Eq. P^ vs. Eq. p^ ]. and the poles in the complex 
k plane arc not altered as compared to the zero field limit 
[cf. Eq. dlSl) vs. Eq. dSO])]. 

We can also intuitively motivate the scaling ^ 9~^ of 
Eq. ([BS]) for small angles. The strong field regime, char- 
acterized by Eq. ([5^ , is compatible with sending the fre- 
quency w to infinity, i.e.. 



2m, 

U! 
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eeB 
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(71) 



while keeping the combined quantity uj sin 9 fn uj9 finite, 
such that still eeB/{uj9)'^ ^ 1. In this limit, the asymp- 
totics (|68)) should also be valid for any reasonable thick- 
ness d of the barrier (cf . Appendix |B]) . The only remain- 
ing dimcnsionless ratios to govern the photon-to-photon 
transition probability are thus given by 



eeB 



eeB 



and 



4m^ 
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(72) 



Taking into account that the leading term of the pho- 
ton polarization tensor in the strong field limit, Eq. (|T8|) . 
is linear in eeB and expecting the transition probability 
to diminish with diminishing mass, we can straightfor- 
wardly argue that the leading (polynomial) contribution 
of Eq. (|68| should at least scale as 9~^: 
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Finally, note that we have derived expressions for 
the photon-to-photon transition probability in all the 
parameter regimes of interest. In particular wc pro- 
vided analytic asymptotics in the limit w 3> 2m and 

d <g; ^ (^) f*^'^ ^^^ zero, weak and strong magnetic 
field limits. In all these cases, the leading contribu- 
tion to the photon transition amplitude turned out to 
depend on the thickness of the barrier d via a global 
phase factor only, rendering the photon-to-photon tran- 
sition probability independent of d. This behavior can 
also be understood intuitively. The typical size of the 
particle-antiparticle loops along the magnetic field lines 
is of the order of the virtual particles' Compton wave- 
length, i.e., ^ m^^ . In the small mass limit, cj 3> 2m 

and d -^ ^ (it) ' ^'^'^ ^'-"" ^ni^'U ^^ ^^ ^'''^ in a situation 
where d is significantly smaller than the Compton wave- 
length of the minicharged particles. Hence, the interval 
of width d which is not available for photons to split or 
recombine into minicharged particles is negligible, and 
the width of the barrier does not affect the asymptotics 
of the transition probability. As described above, beyond 
the threshold for real particle creation in a magnetic field, 
the asymptotics is governed by the conversion of photons 
into real minicharged particles. As these real particles are 
not reconverted into photons, they do not contribute to 
the photon-to-photon transition probability. Of course, 
this process also does not result in any d dependence. 



IV. A LABORATORY EXPERIMENT FOR LSW 
WITH VIRTUAL MINICHARGED PARTICLES 

A. Comparison of LSW scenarios with minicharged 
fermions 

Being the currently most sensitive laboratory LSW 
experiment, the ALPS collaboration already provides 
bounds [13 on minicharged particles through the LSW 
scenario depicted in Fig. [SJ In models where the frac- 
tional charge of the minicharged particles arises from 
the gauge- kinetic mixing of the standard model [/(I) 
with an additional hidden-sector [/(1)h, 'light shining 
through walls' is possible if photons traverse a light- 
blocking barrier through a real hidden-photon interme- 
diate state Hlzl- 

Obviously, the ALPS bounds on minicharged particles 
rely on the existence of hidden photons. Per se, this 
is not unsatisfactory as many models with minicharged 
particles employ a kinetic mixing mechanism [14l llSf . 
such that minicharged particles often come hand in hand 
with hidden photons. On the other hand, the scenario 
as depicted in Fig. [5] (the external magnetic field is not 
mandatory for the process, but amplifies it) does not give 
direct limits on e but rather only on a combined parame- 
ter that involves e together with the minicharged particle 
to hidden-photon coupling, see, e.g., uM- ^y contrast, 
the LSW scenario with virtual minicharged particles, cf. 
Fig. [21 has the potential to provide direct limits on e. For 
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FIG. 5. LSW scenario with hidden photon and minicharged 
particles. A photon (wiggly Une) can oscillate into a hidden 
photon (zig-zag-line) through a minicharged particle loop. Af- 
ter traversing the wall, the hidden photon can be reconverted 
into a photon by the reverse process. From this scenario 
bounds on the coupling e can be deduced for fixed hidden 
photon to minicharged particle coupling. Note that here the 
external magnetic field is not necessary, but favors the con- 
version process. 



this reason it will be instructive to compare the discovery 
potential of our setup to the currently best laboratory 
exclusion limits on minicharged particles 75|, derived 
from PVLAS [23| polarization measurements. These are 



slightly less stringent in comparison to the ALPS LSW 
data, but do not rely on a hidden-photon intermediate 
state. Lastly note that, ALPS adopts = <};(B, k) = 7r/2 
as this both maximizes the transition probability for 
standard processes involving axion-like particles, and fa- 
cilitates the discrimination of scalar and pscudoscalar 
particle species, see, e.g., [l^. By contrast, for mag- 
netically amplified 'tunneling of the 3rd kind' a choice of 
6* « is favorable. 



B. Experimental parameter set and numerical 
implementation 

In the following, we specify several details of a realistic 
setup with rather conservatively chosen parameters in 
order to illustrate the discovery potential of our LSW 
scenario. To be specific, we mainly use parameters of the 
ALPS experiment as a realistic example. For a dedicated 
experiment, we expect that the optics design could easily 
be optimized. 

As in ALPS-I, we consider light of a frequency dou- 
bled standard laser light source, w = 2.33 eV (A = 
532nmj^, which is fed into an optical resonator cavity of 
length L= 8.6m to increase the light power available for 
minicharged particle production. The resonator uses a 
plano-concave design with one plane mirror and a curved 
one with radius of curvature R = 15m. The incident laser 
light is coupled into the resonator via the curved mirror 
and is directed towards the plane mirror, mounted right 



^ Note that for the employment of a regeneration cavity, it can 
be favorable to use the laser in the fundamental mode, i.e., A = 
1064nm (u! = 1.165eV). This possibility is also taken into account 
in Fig. [3 for 6 = 0.001°. 



in front of the barrier. In this way, a stable resonator 
mode is built up in between the two mirrors [43|. For 
Gaussian beams, the beam waist wq in such a resonator 
is given by [76[, 



A 



Wn 



-^L(R-L). 



(74) 



For our parameters, wq ~ 1.12mm. The beam waist 
corresponds to the radius of the laser beam on the surface 
of the plane mirror. Within the resonator, the radius w 
of the laser beam increases as a function of the distance 
D to the plane mirror |7q |. 



w{D) 




(75) 



Thus, whereas the photons are perfectly aligned with the 
optical axis on the surface of the plane mirror, i.e., at 
D ~ 0, the maximum angular deviation from the optical 
axis A6 in the resonator increases with D, 
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The smallest resolvable minicharged particle mass is de- 
termined by several experimental constraints. 

Firstly, it is limited by the extent of the homogeneous 
external magnetic field. More specifically, the spatial and 
temporal homogeneity of the magnetic field should be 
comparable to or larger than the Compton wavelength of 
the minicharged particles. In addition, our setup requires 
that the direction of the laser beam and accordingly its 
cavity can be adjusted such that they are almost perfectly 
aligned with the magnetic field lines. As emphasized 
above, this discriminates our experimental needs from 
those of existing LSW searches such as ALPS, which cm- 
ploys a dipole magnet with 9 = <(B,k) = tt/2. In con- 
sequence, for our purposes, accelerator dipole magnets, 
as used in ALPS, do not seem favorable, as there the ex- 
isting magnet bore extends transversely to the magnetic 
field lines. The ideal magnet in this tunneling experi- 
ment with virtual minicharged particles is presumably a 
solenoid or Helmholtz coil magneo 

As a suitable magnet we have identified a presently 
unused ZEUS compensation solenoid [T^ available at 



^ Notably, next generation heliosopes and direct dark matter 
searches suggest the employment of large scale toroidal magnets, 
cf. 139 . [73 , which offer a large magnetic field region of the or- 
der of several meters. For our LSW scenario, the employment 
of such a magnet would also be possible, particularly as the AT- 
LAS BO test coil allows for a cavity setup along the magnetic field 
lines. This would open up the possibility to probe minicharged 
particles with masses down to C'(10~'^)eV. 
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DESY. It features a bore of 0.28m diameter and 1.20m 
length and provides a field strength oi B = 5T. The field 
points along the bore, and is assumed to be adequately 
aligned on the solenoid's axis (accurate alignment studies 
of magnetic field lines relative to gravity have, e.g., been 
performed in [7a| for a HERA dipole magnet). The field 
strength near the center of the solenoid is expected to be 
sufficiently homogeneous at least over a typical extent of 
the order of the bore's diameter. The wall is installed in 
the center of the bore and the back end of the cavity ex- 
tends into the bore. The angle 9 is adjusted by tilting the 
entire optics assembly relative to the solenoid's axis. As 
discussed below Eq. (fT3|) , the length scale over which the 
field can be considered as approximately homogeneous 
should be comparable to or larger than the Compton 
wavelength of the minicharged particle. Thus, with the 
ZEUS compensation solenoid, access to minicharged par- 
ticle masses down to m > 7 x lO^^eV - corresponding 
to a Compton wavelength of 0.28m - is granted. As the 
minicharged particles interact extremely weakly with the 
magnet material and all the other parts of the setup, the 
diameter of the beam tube and other parameters are not 
relevant. 

Secondly, as discussed above, the smallest resolvable 
minicharged particle mass is also limited by the pair pro- 
duction threshold U!sm9 = 2m. Note, however, that be- 
yond the pair production threshold, instruments sensitive 
to real minicharged particles could possibly be combined 
with our proposal for an increased sensitivity also in this 
region, see, e.g. [80|- In the present study, we do not 
consider such options. 

In practice, the relation for the pair production thresh- 
old directly constrains the smallest resolvable mass via 
the smallest experimentally realizable angle 9. On the 
one hand, the smallest achievable angle depends on the 
magnet bore, as the generation cavity (and possibly also 
the regeneration cavity, see below) needs to fit into the 
magnet bore and thus may not exceed a certain size. On 
the other hand, it must be "large enough" as to allow for a 
precise geometric alignment. Below, we will give results 
for three exemplary angles, namely for 9 = 0.1°, 0.01° 
and 9 ~ 0.001°. Given a precise alignment of the mag- 
netic field lines within a diameter of ©(mm) about the 
solenoid's axis - recall, that wq w 1.12mm -, the un- 
certainty in the adjustment of 9 is expected to be domi- 
nated by the maximum angular deviation from the opti- 
cal axis in the resonator A9. For D = 0.6m, i.e., half the 
length of the ZEUS compensation solid, Eq. ([75)1 yields 
A9 « 0.0007°. Given this small uncertainty in the align- 
ment, angles down to 9 = 0.001° should indeed be exper- 
imentally feasible. In |5l| . we have even overestimated 
this uncertainty. 

In addition, we comment on the thickness d of the 
light-blocking barrier, which is assumed to be perfectly 
refiecting, cf. Sect. |lll For ALPS, the barrier is real- 
ized by a high-quality steel block of width d = 1.8cm = 
9.1 X lO^eV"^ [8l|: As for ALPS the barrier is assumed 
to be traversed by real particles, the photon-to-photon 



transition probability does not depend on the thickness 
of the wall in the limit of vanishing WISP-to-matter cou- 
pling. 

By contrast, our tunneling scenario exhibits an intrin- 
sic sensitivity on the wall thickness: For minicharged 
particle masses corresponding to Compton wavelengths 
smaller than d, the tunneling process is obstructed, 
whereas for Compton wavelengths substantially larger 
than d, the photon-to-photon transition probability even 
becomes independent of d, see also below. To demon- 
strate this, we provide exclusion limits for two different 
thicknesses of the wall, namely d = 1mm = 5.1 x lO'^eV"^ 
and d = Ifim = 5.1eV~ , which are most easily numer- 
ically accessible. However, let us stress that exploring 
large Compton wavelengths, i.e., small minicharged par- 
ticle masses with the considered scenario is experimen- 
tally also possible with thicker barriers, as, e.g., used 
in the ALPS experiment. For minicharged particles of 
larger masses, the experimental challenge can be greater, 
as the minimization of the barrier thickness can become 
an important handle to maximize the transition proba- 
bility. As the discovery potential of our experiment is 
most promising in the small mass regime, no difficulties 
in the "wall" -design should arise. Technically, the nu- 
merical evaluation of the auxiliary function /> becomes 
considerably more involved for thicker barriers, see be- 
low. For this reason, we limit the numerical evaluation 
to the d values stated above. 

The asymptotics derived for the photon-to-photon 
transition probability, Eqs. (|65p and (|67p . are extremely 
helpful here. They do not feature any explicit depen- 
dence on d, and are valid in the parameter regime where 
Lj :$> 2m, and for wall thicknesses d <S^ :^ {^) ■ I* 
is convenient to restate these findings also in a slightly 
different way: For a given laser energy to and wall thick- 
ness d, the asymptotics provide us with the correct re- 
sult of the photon-to-photon transition probability for all 
minicharged particle masses to <g; { ^ , ^/^} ■ As demon- 
strated in Sec. lIVCl below. this will immediately allow us 
to provide results for wall thicknesses of 0(cm), as, e.g., 
employed in ALPS, without the need for any additional 
numerical evaluations. 

We proceed to the evaluation of the observable in the 
proposed LSW scenario, which is the outgoing photon 
rate behind the light-blocking barrier: 



nout = A/" 



R 



7->7 



(77) 



The lever arm for probing small minicharged particle 
couplings is the large number of incoming photons njn 
on the right hand side of Eq. (|77p vs. - in principle 
- single-photon detection for nout- In addition, we in- 
troduce a variable A/", which accounts for the option of 
installing a regeneration cavity on the right hand side of 
the wall. Without a cavity on the regeneration side, one 
has TV = 1. 

Inspired by the experimental possibilities at ALPS, we 
assume an ingoing-to-outgoing photon ratio ni^/ria^t = 
10^^, which already accounts for additional parameters of 
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the experimental realization, i.e., the use of a front-side 
cavity, the effective detector sensitivity and the running 
time, cf. p^ . 

Lastly, we comment on the numerical evaluation of 
Eq. (|77p . Herein a difficulty is the numerical evaluation 
of the auxiliary function /> , contributing to the photon- 
to-photon transition probability P-y-^-y. 

As a function of the integration variable A, the inte- 
grand in /> generically exhibits a highly oscillatory be- 
havior (depending on the parameter set d,uj,m) in all 
relevant limits (with or without external field): Rescaling 

the integration variable as A — > Aw 1 — ^^, the numeri- 
cal evaluation becomes increasingly difficult in the vicin- 
ity of the upper integration boundary A ^ 1 for larger val- 
ues of the oscillation frequency ds/uP — 4™^: For thicker 
walls and masses fulfilling the relation lo ^ 2m, the fre- 
quency of the oscillations is always very large. It becomes 
only slightly smaller for minicharged particle masses in 
the regime w > 2m. In addition also the amplitude of 
the oscillation becomes very large near A ^ 1, especially 
for small m. To circumnavigate this problem, it is help- 
ful to split the A integral at an appropriate cutoff value 
Xc{d,uj,m) which depends on both the set of variables 
used and on the employed numerical integration routine. 
Substituting A — > l/R above the cutoff, we can then 

rewrite the integral as: J^ dX = J^" dX + J^ " dR/R^. 
Thereafter, the numerical routine is capable of treating 
the problematic region. 

In contrast to />, the integrand of /< is never os- 
cillating, but always decaying as a function of the in- 
tegration variable k, cf. Sect. IIIII Nevertheless, for 
the numerical evaluation, a similar rescaling procedure, 

K — >■ kJ ^ 1 for uj < 2m is consistent. For w > 2m, 

the numerical evaluation of the k integral is straightfor- 
ward. 

Let us emphasize, that as we have worked out asymp- 
totics for all relevant limits, we can easily obtain esti- 
mates for different parameter sets, without having to go 
through a full numerical analysis. 



Prospective exclusion bounds for minicharged 
fermions 



The discovery potential of our setup can be quan- 
tified by prospective exclusion bounds for minicharged 
fermions, see Figs. [6]and[7l The colored areas depict the 
parameter regimes that are accessible by the proposed 
experimental setup, as well as different existing bounds, 
see below. The two plots concentrate on the meV mass 
range (Fig. |6]), and on a mass range including very small 
masses (Fig.[7])- 

Fig. [6] displays the mass-coupling plane for 'large' 
minicharged particle masses of m — 10~^cV to m ~ 
2eV. This mass range is obviously irrelevant in terms 
of its discovery potential for minicharged fermions, as all 
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FIG. 6. Exclusion bounds in the fractional-charge mass 
plane for 'large' minicharged particle masses. Without an ex- 
ternal field, the uppermost area, depicted in blue (dark gray), 
can be excluded. The white dashed line is determined with 
Eq. (|50|l . It delimits the excluded area from below. As the 
underlying weak field approximation is limited to ^^ <g 1, 
the corresponding bounds can only hold to the right of the 
yellow dot-dashed line, representing ^^ — 1, where they do 
not significantly improve the zero-field bounds. For any fi- 
nite magnetic field, however, a divergence is found in the 
photon-to-photon transition probability for cu — 2m (cf. main 
text). The area above the orange line, particularly the or- 
ange (light gray) area in the lower left labeled by 'strong', 
can be excluded with Eq. ^^, valid for {^, -^jt^^) > 1, 
when assuming angles 6 < 0.1° and using a wall thickness of 
d — Ifim for easy numerical evaluation. The black dashed line 
defines the area excluded when setting d = 1mm instead. It 
can be seen that the bounds become independent of d below 
m ^ 5 X lO^^eV. Notably, the orange area is also perfectly 
reproduced when employing the analytic asymptotics (|65|) in- 
stead of Eq. (|59p . The green, dash-dotted vertical lines de- 
pict the condition d = i (^)~ (cf. App. |B]) for d = 1mm 
and d = 1.8cm. For a given d, the asymptotics (|65|) become 
trustworthy to the left of the respective line. This theoretical 
prediction is perfectly compatible with the convergence of the 
black and the orange lines to the left of the d = 1mm line. 
Hence, most importantly, all exclusion bounds to the left of 
the d = l.Scm (the wall thickness at ALPS) line are also valid 
for wall thicknesses of up to at least d — l.Scm, even though 
the respective numerical calculations employ d — l/im. As 
a reference, the parameter regime excluded by the PVLAS 
experiment is depicted in purple (gray). 



the regions accessible by magnetic tunneling via virtual 
minicharged particles are already excluded by PVLAS 
[24| polarization data (cf. the lowermost, purple (gray) 
shaded area). We can nevertheless obtain several in- 
teresting physical insights. A numerical evaluation of 
Eq. ([55]) . using d = 1/im, results in the blue colored area 
in Figs. |6] and [T] As already discussed in [52], going to the 
smallest considered masses of m = 8 x 10~'^eV, the low- 
ermost boundary of the excluded area in the absence of 
an external field decreases with a logarithmic dependence 
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and finally reaches a fractional coupling of e ~ 6 x 10~^, 
cf. also Fig. [71 Returning to Fig. [51 the area above 
the white dashed line, is the discovery potential in the 
limit of weak magnetic fields for 6* sa 0, cf. Eq. ([5n|. 
As argued in Sect. IIII El the perturbative calculation is 
only valid as long as ^^ ^ 1 . In Fig. [HI the perturba- 
tive correction can be expected to be trustworthy only to 
the very right of the yellow dot-dashed line, correspond- 
ing to ^^ = 1. Its limit of validity is clearly exceeded 
before it can improve the zero-field bounds. Therefore, 
the perturbative, weak-field result adds hardly anything 
to the new physics discovery potential in the zero-field 
limit already discussed. For the exploration of an hith- 
erto untested minicharged particle parameter regime, the 
strong-field results of Sect. IIII El arc crucial. A distinct 
feature discriminating the zero-field situation from the fi- 
nite magnetic field case is certainly the resonance in the 
photon-to-photon transition probability for w = 1m. For 
our parameters it is found at m = 1.165eV, see Fig. [51 

As argued in Sect. IIII El for ^^ ^ 1, i.e., to the very 
left of the yellow dot-dashed line, the strong field ap- 
proximation, Eq. (j59p . is valid. The additional constraint 
.^ '2 „ 3> 1 is met in the entire parameter space shown 
in Fig. [6l We depict the bounds derived in the strong- 
field limit for two different wall thicknesses: For d = 1/im, 
the area above the orange line, particularly the orange 
(light gray) area denoted by 'strong', can be excluded. 
Correspondingly, for d = 1mm the area above the black 
dashed line can be excluded. For minicharged particle 
masses below m ~ 5 x 10~^eV, the lowermost bound- 
aries of these areas fall on top of each other, and hence 
become independent of the wall thickness d. Going to 
larger masses, the lines d = Imni and d = 1/im are found 
to separate, and the thicker wall always results in worse 
exclusion bounds in comparison to the thinner wall. This 
constitutes a nontrivial, full numerical check of the mani- 
festly d independent asymptotics, derived analytically in 
Eq. (|65|) . For a given d, the asymptotics (|65|) are expected 
to become trustworthy if the condition d <C ^ (it) ^® 
met, cf. also the discussion in Sec. IIVBI Strikingly, this 
implies that the exclusion bounds to the left of the ver- 
tical dash-dotted line labeled with d = 1.8cm, although 
numerically computed for d = 1/im, are also valid for wall 
thicknesses up to at least d = 1.8cm, i.e., the thickness 
of the wall employed in the ALPS experiment. 

Note that the resonance at w = 2m is clearly visible 
in the strong-field results in Fig. [5] It is, however, lo- 
cated in the region where ^^ < 1, such that the strong 
field approximation cannot be expected to hold anymore. 
Rather, the weak field approximation should be valid 
close to the resonance at w = 2m. 

In Fig. [7] we present exclusion bounds in the fractional- 
charge mass plane ranging from minicharged particle 
masses of to = 8 x lO^^eV to m = lO^^eV. The orange- 
colored cone-shaped area depicts the parameter space 
that could be excluded by virtue of Eq. ([5^ for an an- 
gle oi 9 = 0.1° and a wall thickness of up to at least 
d < 1.8cm, cf. the discussion in the context of Fig. [5] 
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FIG. 7. Exclusion bounds in the fractional-charge mass plane 
for 'small' minicharged particle masses, valid for wall thick- 
nesses up to at least d = 1.8cm. All results besides the zero- 
field bounds are derived numerically from Eq. (|59l) . The dark 
blue (darkest gray) area depicts the parameter regime that 
can be excluded in the zero-field limit, cf. also Fig. (6] The 
orange (light gray), cone-shaped area on the very right depicts 
the parameter space that can be excluded when choosing an 
angle of ^ = 0.1°. Thereafter, from right to left, pink (gray) 
and dark red (dark gray) shaded cones correspond to angles 
of e = 0.01° and Q = 0.001°. The lowermost light red (gray) 
shaded cone is also for Q = 0.001°, but in addition accounts 
exemplarily for a regeneration cavity. The centers of the cone- 
shaped areas are marked by dot-dashed lines, denoting the 
respective pair production thresholds, oj sin d = 1m. The 
steep cusps at the pair creation thresholds might be smoothed 
and less pronounced in an actual experimental realization, 
see main text for details. Note that the exclusion bounds 
for smaller angles always include the entire parameter space 
which can be explored for a larger angle. The asymptotics for 
S — 7- is given by the black dashed line, whereas the asymp- 
totics for 0—^0 with a 2nd cavity added is marked by the 
yellow (gray) dashed line. We can also provide for analytic 
asymptotics beyond the pair creation threshold, cf. Eq. (|68|) . 
Such asymptotics are shown exemplarily for Q = 0.1°, as a 
white dashed line. We stress that all these parameter regimes 
fulfill the condition 'i'^^ ' i „ <C 1 for the respective angle 6. 
The dotted, purple line (and area) depicts the limits [75| de- 
rived from PVLAS polarization measurements [23], and the 
blue short-dashed line (and light blue (lightest gray) area) 
refers to the best model-independent cosmological bounds ob- 
tained from CMB data [H. For an angle of 6* = 0.001°, 
and employing a cavity on the regeneration side, magneti- 
cally amplified tunneling with virtual minicharged particles 
can be expected to outmatch these bounds in a mass range 
of 2 X 10"*^ < m < 9 X lO'^'eV. Finally, the pink (lower- 
left) region shows the exclusion bounds for 6 = 0.001° at- 
tainable in a setup with second cavity and infrared laser light 
at oj = 1.165eV, cf. footnote [T] Note that these particular 
bounds are derived from the asymptotics, Eqs. (|65p and 
without evaluating Eq. (|59|) numerically. 



and Sec. IIVBI The center of the cone-shaped area is 
marked by a dot-dashed line, indicating the pair produc- 
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tion threshold uj sin 9 = 2m for the respective incidence 
angle (cf. SeclmEl). 

Following the orange-colored exclusion bounds for 9 = 
0.1°, the cone-shaped areas in pink and dark red (from 
right to left) correspond to the cxplorablc parameter 
space Sit 9 ^ 0.01° and 9 = 0.001°. The lowermost 
light-red cone shaped area corresponds to = 0.001° 
and the additional use of a second cavity. For complete- 
ness, the exclusion bounds for 9 = 0.001° attainable in 
a setup with second cavity and infrared laser light at 
Lo — 1.165eV are represented by the pink (lower- left) re- 
gion. This might constitute a more favorable choice for 
the laser energy in a setting with regeneration cavity, cf. 
also footnote [TJ For simplicity, we have estimated this 
region from the asymptotic formulas (|55|) and (p5)) . Note 
that in this situation also the resonance corresponding to 
9 = 0.001° is shifted according to Eq. ^0^. Remarkably, 
for this choice of parameters, the bounds derived from 
cosmological observations can be overcome in an even 
larger mass range. 

The lowermost boundary of the accessible parameter 
regime in the asymptotic limit — >■ is depicted by 
the black dashed line. The yellow dashed line shows 
the same asymptotics, but in addition assumes resonant 
regeneration with a second cavity installed behind the 
wall, employing Af = 10^. As the transition probability 

(strong) 



P-)!];™"^^ ~ e^, cf. Eq. dSni), employing the second cav- 
ity yields almost an order of magnitude better exclusion 
bounds. Note that in order to achieve the same effect 
by tuning the magnetic field strength, one would corre- 
spondingly need to enhance the magnetic field by a factor 
of a/IO^ as it enters Eq. ([55]) by its square. This would 
demand for magnetic fields oi B ^ 1.5 x lO'^T, corre- 
sponding to field strengths that can only be obtained in 
highly focused lasers in a laboratory. However, there of 
course the extent of the field itself would be limited to 
a smaller region and requires to go beyond the homoge- 
neous field approximation used here. Thus, a large effort 
on the optics side seems more beneficial than requiring 
increasingly stronger magnetic fields. 

To allow for a comparison of our results with current 
experimental limits. Fig. [7] includes the exclusion limits 
75 1 as provided by PVLAS polarization measurements 
2J] (upper dotted purple line and area). The light blue 
area, limited from below by the blue dotted line, de- 
picts the best model- independent cosmological bounds 
obtained from CMB data |42| . 

To summarize: even with conservatively chosen exper- 
imental parameters, we find that the LSW scenario with 
virtual minicharged fermionic particle-antiparticle states 
could improve the limits provided by the PVLAS col- 
laboration below m < 10~^eV, and even outmatch cos- 
mological bounds derived from CMB data below m < 
5 X 10"^eV for extremely small angles — > 0. Em- 
ploying a cavity on the regeneration side, these values 
can be even improved to outmatch PVLAS bounds for 
m < 2 X 10~'^eV and to overcome cosmological bounds 
at m < 9 X 10-^eV. 



V. CONCLUSIONS 

In this paper, we have generalized the recently pro- 
posed light-shining-through-walls scenario via virtual 
minicharged particles [52[ to account for a homogeneous 
external magnetic field. Moreover, we have demonstrated 
how this particular LSW scenario can be employed to 
constrain the parameter space for minicharged fermions 
in a dedicated laboratory experiment. Whereas in the 
absence of an external field, the tunneling scenario with 
virtual minicharged particles has only little sensitivity 
to a new-physics parameter space J52l |. we have shown 
that an external magnetic field can significantly enhance 
its discovery potential. Notably, magnetically amplified 
'tunneling of the 3rd kind' even allows, in a rather small- 
scale laboratory setup, to explore a parameter regime 
which has been inaccessible so far. 

To this end, we have suggested a simple experimen- 
tal setting, involving only present-day technology. In 
contrast to standard LSW setups, typically employing 
9 = <(B,k) = 7r/2, our experiment requires the probe 
photons to enter at a preferably very small angle 9 with 
respect to the homogeneous magnetic field lines. To- 
gether with a resonant regeneration cavity, it has the po- 
tential to outmatch the currently best direct laboratory 
bounds on minicharged particles in the sub-meV mass 
range, as provided by the PVLAS collaboration. In the 
range of minicharged particle masses below O(10~'*)cV, 
it could explore - and even overcome - a parameter 
regime so far only accessible in large scale cosmological 
observations. 

We note that similar ideas could also be considered on 
terrestrial or astrophysical scales, as magnetic fields of 
larger extent might provide access to extended regions of 
the minicharged particle parameter space. Also, a modi- 
fication of this setup beyond the optical regime might be 
favorable in the search of higher-mass minicharged par- 
ticle dark matter [j, |j] . 

From a formal point of view, our analysis substan- 
tially builds on analytic insights into the QED photon 
polarization tensor at one loop accuracy, while retaining 
its full momentum dependence. Within standard LSW 
scenarios, the barrier is traversed by means of real par- 
ticles, whose energy and momentum is fixed by energy- 
momentum conservation. By contrast, the tunneling sce- 
nario via virtual particles genuinely requires to account 
for the entire range of allowed four momenta within the 
particle-antiparticle loop. This significantly complicates 
the corresponding calculation, but also gives rise to sev- 
eral novel features, as compared to conventional LSW 
scenarios. 

We have derived analytical expressions for the photon- 
to-photon transition probability in the limiting cases of 
both, weak and strong magnetic fields, corresponding to 
{ll^^lf } « 1 and i^,^;^} » 1, respectively. 
An exemplary numerical evaluation of these expressions, 
using experimental parameters which can be feasibly im- 
plemented, has allowed us to provide exclusion bounds 
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for fcrmionic minichargcd particles in a wide range of 
the fractional-charge mass plane. 

In addition, we have analytically extracted simple 
asymptotics for the photon-to-photon transition proba- 
bility in the most relevant limiting cases, such that the 
obtained results can be straightforwardly adopted to dif- 
ferent experimental parameter sets. 

In summary, the relatively small experimental effort 
required to realize the "Tunneling of Light via Loops" 
experiment proposed by us compares favorably with the 
prospect of exploring a parameter space for minicharged 
particles hitherto inaccessible in the laboratory. An im- 
plementation of this setup thus seems highly worthwhile 
and will put the existence of minicharged particles fur- 
ther to the test. 



leading contribution is ~ eeB [cf. Eq. ((T8)) ]. the photon 
polarization tensor for scalar minicharged particles in the 
strong-field limit is suppressed and starts to contribute 
only at subleading order |66| . The reason for this qualita- 
tive difference lies in the fact that the bosonic fluctuation 
spectrum in the presence of a magnetic field does not ex- 
hibit a near-zero-mode. We conclude that a Pauli-type 
spin-field coupling is essential for the physics discussed 
in the main part of this work. Hence, given the same ex- 
perimental parameters, we expect the exclusion bounds 
on scalar minicharged particles, achievable via 'tunneling 
of the 3rd kind' to be much weaker. 



Appendix B: The regime of validity of the analytic 
asymptotics in the small-mass limit 
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Appendix A: Scalar minicharged particles 

A perturbative expansion of the photon polarization 
tensor for scalar minicharged particles [S^] analogous to 
Eq. dH]) yields 



nf 
nr 



!!■ 



sc(0) 



jjsc(2) 



+ 0{ieeBr). (Al) 



The zero-field expression reads 
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and the contribution at O ((ee_B)^) is given by 
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with $0 defined in Eq. ([T5| . Implementing the 
limit eeB — >■ oo in the polarization tensor for scalar 
minicharged particles, the first nonvanishing contribu- 
tions are at most logarithmic in eeB. Hence, in com- 
parison to the photon polarization tensor for fermionic 
minicharged particles in the strong-field limit, whose 



In this appendix, we want to comment on the range of 
validity of the various analytic asymptotics in the small- 
mass limit, Lo ^ 2771, derived in the main text. The 
asymptotics turned out to be independent of the thick- 
ness d of the barrier. Here we study the d regime where 
these asymptotics can be considered as trustworthy. 

To derive the asymptotics, we basically expanded the 
integrand in /> around A = 1, and performed the A in- 
tegration thereafter. We thereby also approximated in 
particular the exponential part in the integrand as 



cdX 



(Bl) 



The latter term, being a pure phase factor, canceled 
in the determination of the photon-to-photon transition 
probability, as it involves an overall modulus-squared, cf. 
Eq. ([55)1 . Obviously the approximation in Eq. (|B1[) is 
well-justified for Lod ^ 1. as to leading order in ojd <C 1, 
both exponentials in Eq. (|Bip are equal to one. Its regime 
of validity in the opposite limit Lod 3> 1 is less clear. 
Hence, we subsequently focus on the limit tod ^1. In 
order to estimate when this approximation is justified, 
we study the following integral. 
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for 77 = 3,5,7,..., both by setting A = 1 in the expo- 
nential from the outset, and by explicitly performing the 
integration over A. Note, that Eq. (jB2[) exhibits the ba- 
sic features of the integrals encountered in deriving the 
asymptotics, and can be integrated explicitly. Setting 
A = 1 in the exponential, we obtain 
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whereas in the other case, we first reformulate the cor- 
responding expression with the help of the substitution 
A ^- 1 — /. This results is 
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We now focus on the evaluation of the integral in 
Eq. (|B4[) . For n integer and odd, we obtain 
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with integration constant C. The incomplete gamma 
function r(a, x) has an exact series representation (for- 
mula 8.354.2 of HI), 

r(a,x)=r(a)-^ \,/^/^,^, , a^-No, (B6) 



A:=0 



kl{a + k) 



and its asymptotic expansion for |x| — > cxj is given by 
(formula 8.357 of H) 
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Hence, Eq. (jB4p can be cast in the following form 
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The fc = term in the sum in Eq. (jBlOp agrees with the 
last line of Eq. (|B3p . For n = 3, 5, 7, . . . , and given that 



dm < 1 <^ d<— , (Bll) 

LJ m \ LU J 



it is the leading contribution to Eq. (|B10|) . If Eq. (|B11|) 
is met, both terms due to k > 1 in the sum, and the 
term in the second line of Eq. (jBlOp . become negligible 
as compared to the fc = term. In this limit, Eqs. (jB3[) 
and (jBlOp indeed coincide. Hence, the various asymp- 
totics derived for cu ^ 2m in the main text should be 
trustworthy for thicknesses d of the barrier fulfilling the 
condition rf < — (— )~\ 



with the additional requirement — ^ < argx < ^. 
Equation (jB6p implies 



I I e "'"''' n 



cud /"/2 



-f^(^c.d)"/^[^(-t)-^(-t,^a;dO] 



ln/2 Z^ 



{-ifiudlf 



(B8) 



K2^t^^fc!(fc + l-t) 

We use this expression to evaluate the integral in Eq. (jB4[) 
at the lower integration limit. As we are interested in 
the limit ujd 3> 1 here, at the upper integration limit we 
rather employ Eq. (jB7p . yielding 
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Appendix C: /> ™"^ beyond the pair creation 
threshold 



As pointed out in Sec. IIII C| the pole at A = cos 9 — 
irj becomes relevant in any explicit evaluation of the A 
contribution, if the condition in Eq. p4p is met. It is 
then convenient to rewrite the integral in Eq. ([56]) in a 
manifestly rj independent way. This can be achieved, e.g., 
with the help of an integration by parts. Thereafter, we 
obtain the following rather lengthy, but rj independent 
and finite expression, 
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where we introduced 
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Equation (jCip is not only essential for the purpose of 
explicit numerical calculations in the regime (j34p . it also 
allows for the extraction of an analytical asymptotics for 
the photon-to-photon transition probability in the limit 
w > 2m and d <^ ^ (ir)~^ (^f- App.[B|). 

An expansion in — <C 1 of the terms in Eq. (|C1|) 
featuring no integrations over the parameter A results in 
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We now consider the contributions to Eq. (|Cip involving 
integrations over A. The first integral yields 
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and (cf. the discussion in the main text), 
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Moreover, we obtain 
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and 



C5 = COS ( 



= COS! 



dA 



[g(cos6>)-g(A)] ln|cose-A| 
2 \ 3/2 



4m^ 



dA 



-g(A) ln|cose-A| 
, „ 4m^ 



o 



"2V2(l-|-cost 



4m' 



2\ 3/2 



1 A ln|cos Q — X\ 
^^ (l-A)3/2 



{m") 
<^((^)°) 



e^co.«cos6^1n(l-cos0)^^^^^^^oy ^^^^ 



(1 + cos 61) sin'' 9 



as well as 
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Hence, the leading behavior of Eq. (IC1[) in the hmit 
^^ ^ 1 is given by 
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The terms ^ (^^) ^, encountered in some of the con- 
tributions Cj cancel, resulting in an overall logarithmic 
dependence on the ratio ^^. 
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